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Introduction 

Finite-difference techniques are applied to solve 
field problems on grids which are ordered sets of dis- 
crete points in a coordinate system. Solution accu- 
racy and computational expediency depend on the 
grid as well as the finite-difference technique. For in- 
stance, the accurate representation of boundary con- 
ditions and the resolution of solutions on regions of 
rapid change such as boundary layers, shocks, and 
separations in flow fields require a grid with partic- 
ular characteristics. First, the grid should adapt to 
the physical boundaries so that boundary conditions 
can be readily applied. Second, the grid should be 
concentrated in regions of rapid change to accurately 
compute the solution there. Third, from a computa- 
tional point of view it is desirable that the grid be 
uniform, the boundaries enclose a rectangle in two 
dimensions or a rectangular parallelepiped in three 
dimensions (fig. 1 ), and the exterior boundaries cor- 
respond to physical boundaries. If this is possible, 
overall computer program logic for the application 
of a finite-difference solution algorithm can be mini- 
mized and the process can be kept highly repetitive 
over the entire grid. 

A direct (algebraic) approach to grid generation, 
for which an explicit functional relationship between 
the computational domain and the physical domain 
is known, has the advantage that changes to the grid 
are direct and are rapidly obtained. This report de- 
scribes such an approach. It also describes the use of 
an interactive computer program for constructing a 
direct functional relationship between the computa- 
tional domain and the arbitrary, simply connected, 
two-dimensional physical domains such that bound- 
aries in the two domains map into each other. The 
technique is called the “two-boundary technique” 
and is described in general in references 1 to 4. The 
outlining features of the technique are derivatives of 
transfinite interpolation (ref. 3) and the technique 
can also be, in part, derived from the general mul- 
tisurface equation (ref. 4). The two-boundary tech- 
nique is best applied in an interactive environment. 
An interactive program based on the technique has 
been written, and its usage is described herein. 

Background 

For a large class of problems, an approach to grid 
generation which tends to satisfy both the accuracy 
and expediency requirements is to transform the gov- 
erning equations from the original defining coordi- 
nate system (referred to as the physical-coordinate 
system) to a rectangular computational-coordinate 
system. The particular region of a coordinate sys- 
tem where a grid is defined is referred to as a domain. 


The computational domain can be expressed by 

0 < £ < 1 
0 < 77 < I 
0<f < 1 

where £, 77 , and f are computational coordinates. (A 
list of symbols and abbreviations used in this paper 
appears after the references.) A transformation from 
the computational domain to the physical domain 
can be expressed as 

X = a ;(£,»?,?) ' 

y = y{t,v,$) > (1) 

where x ) y, and z are physical coordinates. Simi- 
larly, a grid in one domain can be mapped into a 
grid in the other domain (fig. 2 ). When the transfor- 
mation maps boundaries in the computational do- 
main into boundaries in the physical domain, the 
term “boundary-fitted coordinate system” (ref. 1 ) is 
used to describe the transformation. 

The transformation of the governing equations 
implies that derivatives with respect to the physi- 
cal coordinates must be transformed to derivatives 
with respect to the computational coordinates. Ref- 
erence 1 describes the transformation of the govern- 
ing equations of fluid flow for which the chain rule of 
calculus is applied. The basic result of the transfor- 
mation is inclusion of the Jacobian matrix 



in the transformed equations of motion. Thus, for the 
application of finite-difference techniques, the pri- 
mary function of grid generation is the determination 
of grid derivatives at grid node points. Normally, the 
process is to express the physical grid as a function 
of the computational grid and then to differentiate 
to obtain the inverse Jacobian matrix, which is itself 
inverted to produce the Jacobian matrix (ref. 1 ). 

The application of the two-boundary technique 
can be thought of as two distinct processes. The first 
process is the definition of boundaries (referred to as 
the bottom, top, left, and right) and the functional 
description of the inscribed physical domain. The 
second process is the determination of the properly 
spaced grid in the physical domain which corresponds 
to a uniform grid in the computational domain. 

A convenient way to represent the boundaries 
is to specify subsets of points and an interpolation 



procedure to specify arbitrary boundary points. The 
choice of a boundary interpolation procedure can be 
dependent upon many things. For the interactive 
grid-generation program, it is assumed that the sub- 
set boundary definitions are densely specified and lin- 
ear interpolation is sufficient for arbitrary boundary 
definition. Accuracy of the boundary representation 
is directly controlled by the number and location 
of the subsets of points. The essence of the grid- 
generation technique is connecting distributions of 
points on the bottom and top boundaries. Here sim- 
plicity and versatility are achieved by use of Hermite 
cubic polynomials as connecting functions. If only 
endpoints of the cubics are specified, the connect- 
ing functions are straight lines; however, if the po- 
sition and the derivatives orthogonal to the bottom 
and top boundaries at the endpoints of the cubics 
are specified, an effective curved connecting function 
is defined (fig. 3). A distribution of points along 
the connecting function forms the initial grid. In 
the event that side boundaries to the domain are 
specified, the connecting curves between the bottom 
and top boundaries must be made to conform. Lin- 
ear blending functions are applied to the initial grid 
to create a final grid which conforms to the side 
boundaries. 

The second process, and the one most depen- 
dent upon interactive graphics, is the distribution of 
grid points along the bottom and top boundaries, 
the distribution of the magnitude of normal deriva- 
tives at the boundaries, and the distribution of grid 
points along the connecting function. A smooth- 
cubic-spline technique (ref. 5) has been adopted to 
control these distributions. The degree of smooth- 
ness is a function of user-specified parameters. The 
reason for using a smoothing process is that the ver- 
satility and continuity characteristics of cubic splines 
can be employed and the oscillation problem that can 
occur with unsmoothed cubic splines can be avoided. 
The control of the distribution of grid points is based 
on establishing a parametric variable such as approx- 
imated arc length along a curve (i.e., a boundary 
curve or a connecting curve), normalizing the para- 
metric variable with respect to its maximum value, 
expressing each of the coordinates of the curve as 
a function of the parametric variable, and express- 
ing the parametric variable as a function of a com- 
putational coordinate which is defined on the unit 
interval (i.e., 0 < £ < 1). The unit square is a 
control space (fig. 4) for the distributions along the 
curves. A discrete uniform distribution of the com- 
putational coordinate maps into an arbitrary distri- 
bution of the physical coordinates. Using interactive 
graphics to digitize a few points on the unit square for 
each control function allows for arbitrary control of 


the physical-grid distribution. For the two-boundary 
technique applied in two dimensions there are seven 
distributions to be determined. The primary pur- 
pose of the interactive program described herein is 
the determination of these control functions. 

The Two-Boundary Technique 

The two-boundary technique for grid generation 
is a methodology for establishing the mathematical 
expression (eqs. (1)) relating the computational do- 
main to a physical domain. The methodology sepa- 
rates the boundary-grid definition from the interior- 
grid definition. The boundary grid is first defined, 
and then the interior grid is defined as a function 
of boundary position, boundary derivatives, and an 
independent variable t (0 < t < 1) spanning the 
two boundaries. Because of this approach, the alge- 
braic expression of equations (1) in two dimensions 
is rewritten as 

x = x(X B (r),Y B (r), X T ( S ),Y T (s), t)\ 

V = v{X B (r),Y B {r), X T (s),Y T (s ), t) / U 


where 0 < r < 1, 0 < s < 1, 0 < i < 1, 
X#(r), Yg(r) are coordinates on the bottom bound- 
ary, and Xj’(s), Yj-(s) are coordinates on the top 
boundary. The parametric variables r, s, and t 
are expressed as functions of the computational 
coordinates. 

The boundaries are described by ordered sets of 
points 



where Nq and Nx are the numbers of defining points 
in each set. They are parameterized by computing 
the approximate arc lengths from the first points 
through all the points in the boundary sets. That 
is, 


f l = f- 1 + [{X { B - X^ 1 ) 2 

, 1/2 


+ (Yb - y B l ) ] 
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The approximate arc length corresponding to each 
point is normalized with respect to the maximum 
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value for each boundary by 
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With the sets $1, $2> and $4, arbitrary points 
and derivatives along the boundaries are described 
by 


Xq{v), Yb( v ) interpolation of r into $1, $2 

Xj’(s), Y t {s) interpolation of s into $3, $4 


derivative interpolation of r 
into $1, $2 


dX T {s) dY T (s) 
ds 5 ds 


derivative interpolation of s 
into $3, $4 


For the description herein all boundary interpola- 
tions are linear; however, in general the interpolation 
procedure along the boundaries is arbitrary. 

The cubic function (refs. 1 and 2) connecting the 
two boundaries is 


* = X B (r) on(t) + X T (s) 02(f) + 03(f) 

dr 

+ Q(s)~^- a 4 (t) 

ds 

V = Y B (r)cn(t) + y T (s)a 2 (f) - p(r)^YL a 3 (t) 

dr 

-Q( S )^a 4 (t) 


(3) 


where 


ai(£) — 2t 3 — 3t 2 + 1 
^2(0 “ — 2t^ -f- 3£ 2 
a s{t) — — 2t 2 -\- 1 

0:4 (t) = £ 3 — t 2 


(0 < t < 1) 


P(r ) is the magnitude of the normal derivative along 
the bottom boundary, and Q(s) is the magnitude of 
the normal derivative along the top boundary. 

The cubic curve in equations (3) is specified by: 
(1) the coordinates on the bottom and top bound- 
aries which are functions of approximate-arc-length 


variables, (2) the derivatives of the bottom and 
top coordinates with respect to the approximate- 
arc-length variables, and (3) the magnitudes of the 
derivatives at the bottom and top points. The blend- 
ing functions ai(i), <22 W, <*3(0» an d 0:4^) are shown 
in figure 5. 

The magnitudes of the derivatives affect the or- 
thogonality of the cubics relative to the boundary 
curves and are expressed as functions of the approxi- 
mate arc lengths. Increasing the magnitudes extends 
the orthogonality of the cubic away from the bound- 
ary curves, and excessively large values of the mag- 
nitudes can cause the cubics to intersect themselves 
between the boundaries, which is unacceptable for a 
grid. 

The approximate arc lengths along the bottom 
and top boundaries are expressed as a function of 
the computational coordinate £. There are four 
distribution functions that determine the spacing of 
the cubics. They are 


r = fl (0 ' 

s = 1 ( 4 ) 

P(r) = K B f 3 (0 I W 

Q{s) = K t } 4 ( 0 . 

The functions /i(£), /2(f), /3(f). and MO are 
called control functions, and their determination is 
the second phase of the grid-generation process. The 
constants Kg and Kj> scale the magnitudes of the 
derivatives at the bottom and top boundaries. Fig- 
ure 6 depicts the process for determining arbitrary 
boundary points and derivatives from the initial tab- 
ular descriptions and the connecting function. 

A distribution of the variable t between 0 and 1 
specifies a distribution of points along the curve for 
which the first and last points are at the top and bot- 
tom boundaries (fig. 3). The variable t could be made 
a direct function of the computational coordinate rj. 
but instead it is made a function of the normalized 
approximate arc length along each cubic and then 
the approximate arc length is made a function of r). 
This step requires that equations (3) be used twice to 
compute a grid and is included because the applica- 
tion of side boundaries, to be discussed subsequently, 
is dependent on an arc-length parameterization. The 
variable t is made an empirical function of normalized 
approximate arc length by inversely computing the 
approximate arc length as a function of t at uniform 
steps in equations (3). Values of t versus normalized 
approximate arc length are determined from 

'i = j^3T tf=l,2 Np) 


3 



Uj 


-l + {(Xj -tf-x) + (yj - Vj- 1) 1/2 
[u\ — 0 ;j — 2,3, 


and 


u j — 




(0 < < 1 ) 


where Xj and yj are computed from equations (3). 
For each cubic there is a set 


# — {Uj,tj} 3 j _ i p 

The number of points Np along the cubic connecting 
function is chosen to yield a sufficiently accurate 
approximation of the arc length. Since u and t are 
monotonically increasing functions of each other, the 
values of ^ can be interpolated for a distribution 
of the variable t corresponding to a distribution 
of the approximate arc length. The normalized 
approximate arc length is expressed as a function of 
the computational coordinate 77 by u — 75(77). 

Given the boundary data, there are now five 
control functions that are necessary to map a uniform 
computational grid into a physical grid by use of the 
two-boundary technique. Thus, for 


^ 7 ) = inT (J r = 


is computed. Figure 7 depicts the process of com- 
puting a grid with equations (3) and the boundary 
data. j-m 

At this point, {x(l, J), y(l, J)} and {x(N, J), 

y(7V, define the left- and right-boundary 

grid points. However, prespecified left- and right- 
boundary curves are defined by the sets 

< X k v k X k=N L 

\*L> y Lf k =l 
i X Ri Y R } k=1 R 

where Np and Np are the numbers of defining points 
in each set. The side boundaries are parameterized 
by computation of the approximate arc lengths from 
the first points through all the points in the boundary 
sets. That is, 

v k = v k ~ 1 + [(x£ -X k L - 1 ) 2 

+ (yfc-K*- 1 ) 2 ] 172 (v 1 =0;fc = 2,3 ,...,N l ) 

w k = w k ~ i + [(X& - X £ -1 ) 2 

+ (Yr ~ Yt 1 ) 2 ] V2 (& 1 = 0; k = 2, 3, .... N r ) 


??(J) = ^T (J = 1,2, 

which is a uniform distribution of the computational 
domain, where N is the number of grid points in the £ 
direction and M is the number of grid points in the 77 
direction. Distributions of normalized approximate- 
arc-length parametric variables, 
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1=1 


{*/} 


I=N 

1=1 


i U j ) 


J=M 
J= 1 


are computed with the functions /i(£), /2(C)? and 
75(77), and the distribution 

J=M 

1=1 


is computed with the functions /3(£), /4(£), and 
equations (3). Finally, again through use of equa- 
tions (3), the grid 




J=M 

I=N 


J= 1 
1=1 


The approximate arc length corresponding to each 
point is normalized with respect to the maximum 
value for each boundary by 
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Grid coordinates for the left and right boundaries are 
obtained by interpolating u(J) = fa{r}(J)) into the 
sets 0i, 0 2 , 03, and ©4 such that 


%l{J) = Interpolation of u(J) into 0i = {v k t 
Vl{J) = Interpolation of u(J) into 02 = 
xr(J) = Interpolation of u(J ) into 03 = {w k y 
Vr(J) = Interpolation of u(J) into 04 = {w k > 


4 



The final grid that conforms to the specified left 
and right boundaries is obtained by blending, where 


*(/, J) = *(/, j) + [x L (J) - *(i, J)]{1 - / e [€(7)l} ' 
y(I,J) = »(/, J) + [y L {J) - 9 (1,J)]{1 - feltm * 

+ lvn(J)-y(N,J)}{f 7 [m} 


(5) 


where I = 1,2, ..., N and J = 1, 2, ..., M. The con- 
trol functions /e[£(/)] and /7[£C0] specify the blend- 
ing of the side boundaries with the grid that is ob- 
tained through use of the bottom and top boundaries 
(fig. 8). 


into the control functions. This problem is avoided 
by using a smoothing cubic-spline technique and 
specifying the amount of smoothing as well as the 
digitized points in the unit square. The technique 
that is used is described in detail in reference 6 and is 
described in less detail in reference 5. The features of 
the application of smooth cubic splines are outlined 
below. 

For grid spacing control, the independent variable 
0(0 < 0 < 1) represents a computational coordinate 
and the variable 0(0 < 0 < 1) represents a control 
function. A discrete set of points 

•4 = {fc<Cf° 


Grid Spacing Control 

The second phase of the algebraic grid-generation 
technique is the control of the spacing of grid points. 
For uniform distributions of f and rj , the functions 

MO* MO* MO* MO* MO* MO* and MO 

along with the constants Kq and K? determine 
distributions of the parametric variables r, s, and u 
and the derivative magnitudes P(r) and Q(s). The 
functions used to control the spacing of grid points 
can be analytical. For example, the function 

Mv) = 

— 1 

would concentrate grid points close to the bottom or 
the top boundary depending on the magnitude and 
the sign of the constant K. 

For a general grid- generation procedure, greater 
flexibility in the control of grid spacing is needed than 
can be obtained from specific analytical formulas. A 
general approach is the application of cubic-spline 
functions constrained to be inside the unit square. 
The essence of the control domain (unit square, see 
fig. 4) is that the abscissa corresponds to the per- 
centage of grid points and the ordinate corresponds 
to a particular control function which, in turn, re- 
lates to the geometric definition of the physical do- 
main. A control function can be specified by digitiz- 
ing a few points in the unit square and then applying 
the spline continuity conditions. The functions /i(£)> 

MO* MO* MO* and MO P ass through the origin 
and monotonically increase to the point (1,1). The 
functions MO and MO mus t be single valued but 
otherwise are free to be anywhere in the unit inter- 
val. The derivative magnitudes P(r) and Q(r) are 
obtained by multiplying the constants Kq and K j- 
by MO and MO, respectively. 

The difficulty with using cubic splines for control 
is that oscillations can be inadvertently introduced 


specifies a cubic spline. The conditions on the points 
for the functions MO, f 2 {0, Mv), MO, and f 7 (0 
are 

0 < 02 < 03 < ... < 1 
0 < 02 < 03 < - < 1 

For the functions /3(f) and MO the conditions are 
0 < 02 < 03 < ••• < 1 
0 < 02 < 03 < - < 1 

A cubic spline in the (0,0) plane with the knots 

{<t>iA)Zi C is 

m = {M0i)}£f c_1 

= Wi + M0 - 0t) + C*(0 - 0i) 2 
+ di{ 4 > - 0*) 3 }]~f c 1 


where 


V'iifa+l) — /4i+lW>i+l) 

Ark ~ , , (^«+l) 


d<t > 2 ^ ,+1 ^ 


d 2 n 

~~d<P 


i+1 


(01+ 1) 


(i = 1, 2, ..., Nc~2\ 0 < <j>i < 1) 


The coefficients 1 C 1 are undeter- 

mined parameters whose solutions define F(<j>). The 
objective is to find the coefficients which minimize 
the integral of the second derivative squared 

/V"(0)] 2 d0 

Jo 

subject to the constraint 


N C 

£ 


2=1 


m)-o l 


< C 


5 



where C is a positive constant specifying the extent 
of smoothing and 86 i is the allowed deviation of the 
spline function from the ith ordinate The restated 
objective is to find the smoothest cubic spline passing 
within the bounds 

6i — 86{ < 9{ < 6{ + 89 j 

b = 

69\ = 89jy c = 0 

where 89 ^ represents the maximum derivations of the 
spline function from the specified ordinates 6^ (fig. 9). 

The method of Lagrange multipliers from the 
calculus of variations is used to find the parameters 

The solution algorithm for spline smoothing can be 
found in reference 6. The techniques exist in subrou- 
tine form, where the sets A and B and the constant 
C are input. 

The set A is obtained by digitizing points on the 
unit plane, the constant C is set equal to Nq — 1, 
and the set B is obtained from the linear relation 

66i = m 0 + mi f ^ ,+1 ~ 6 *~ l \ (i = 2,3,...,iVc - 1) (6) 
\<Pi + 1 

Using this relationship allows rapid choice of the de- 
viations which prescribe the amount of smoothing. 
Choosing tuq and mi to be zero results in B = 0, 
and the control function F(<p) is a cubic spline fit to 
the set A. Cubic splines are subject to oscillations 
which are easily observed by plotting the derivative 
functions. Choosing mo and mi to be large results in 
large allowable deviations in the set B, and the con- 
trol function F(<j>) is a straight line. The objective is 
to allow just enough dispersion in set B to achieve 
a smooth control function which satisfies the condi- 
tions described above and produces the desired grid 
concentrations. This process is ideally suited for in- 
teractive computer graphics, where the set A can be 
digitized with a cursor, the constants mo and mi are 
alphanumeric input, and the results can be rapidly 
evaluated. 

Interactive Algebraic Grid Generation 

A computer program based on the algebraic tech- 
nique described above is computationally simple be- 
cause only boundary coordinates and control points 
on the unit square are required input. Linear in- 
terpolation, the connecting function, and the spline 
smoothing technique are the mathematical compo- 
nents. An interactive program, however, is more 


complex than a noninteractive program because the 
user and the program must communicate through 
questions, responses, and graphical displays. Also, 
fault tolerances must be coded so that interactive in- 
put errors will not cause the program to abort. 

The program described herein and called TBGG 
(two-boundary grid generation) is coded in FOR- 
TRAN V and, in its present form, runs on the 
Prime 750 computer and the Control Data CY- 
BER 170 series computers. There are two versions 
of the program — one for the Precision Visuals, Inc., 
DI-3000 graphics library (ref. 7) and one for the 
Tektronix, Inc., PLOT-10 graphics library (ref. 8). 
The DI-3000 version requires 231 000 words of mem- 
ory and the PLOT- 10 version requires 170 600 words 
of memory. With minor modifications, the program 
should run on other computers similar to the ones 
mentioned above. This section describes the program 
input and interactive use. 

There are two forms of noninteractive input to 
the program. In the first form, the boundary 
data are presented to the program in an input file 
(DATANEW, see table I) and a grid based on linear 
control functions and no orthogonality is displayed 
with options for interactive development of the grid 
spacing. The second form of data input is a restart 
file (RESTART) created at the end of a previous in- 
teractive session. The restart file allows the continu- 
ation of the interactive development of the grid spac- 
ing. The output from the program is a free- formatted 
file called GRIDOUT, which contains the grid coor- 
dinates, and an unformatted file called RESTART, 
which contains the parameters and data for the last 
interactively developed grid. 

Program Usage 

Before operating program TBGG to generate a 
grid, a file called DATANEW, which contains the 
problem title, the number of points on each bound- 
ary, and the x- and ^-coordinates of each boundary 
point, must be created. The file DATANEW is read 
by the program when “NEW” followed by a carriage 
return (CR) is the interactive response at the initi- 
ation of the program. If “OLD” followed by a CR 
is the interactive response, the file RESTART (gen- 
erated in a previous session) is read. Examples and 
their associated DATANEW files are described in the 
appendix. 

The output file RESTART is an unformatted 
file and contains the grid spacing information devel- 
oped during an interactive session. The objective of 
RESTART is to recover exactly the conditions that 
existed in a previous interactive session. The file can- 
not be edited, and the specific format and order of 
the data are not described herein. 



After the DATANEW file or the RESTART file 
has been read, a grid and a set of options are pre- 
sented on the terminal screen. The nature of the 
first grid display depends on whether DATANEW or 
RESTART is used for input. If DATANEW is used, 
uniform distributions are assumed for /i(£), /2(C)) 
75(77), /e(£), and /7(f). The orthogonality constants 
K q and K t are set to zero. If side boundaries are not 
specified or if they are straight lines, the curves con- 
necting the bottom and top boundaries are straight 
lines and the spacing between grid points is uniform. 
If RESTART is used, the grid is based on the values 

/i(0, / 2 (C)) fsiOi / 4 (C)) h(v), h(0, M£). K b , 

and Kx previously developed. 


Positive values of the parameter K produce relatively 
small first derivative magnitudes of the control func- 
tion near the point (0,0), and negative values of the 
parameter K produce relatively small first derivative 
magnitudes of the control function near the point 
(1,1) (fig. 12). The variables r and f are used lo- 
cally to represent the independent and dependent 
variables for either of the exponential options. For 
the purposes of the program the range of the param- 
eter K is 0.001 < K < 50. 

Local option 3 uses a combination of two expo- 
nentials to generate nominal control points on the 
unit square. The equations are 


General Options 

The general options are displayed with the grid 
(fig. 10), and they are briefly described in table II. 
An option is invoked by typing the option number 
followed by a CR. 

Boundary and Connecting Function Control 
Distributions (General Options 1, 2, 3, and 4) 

General options 1, 2, 3, and 4 are used to spec- 
ify the distribution of grid points along the bottom 
and top boundaries, the concentration effect rela- 
tive to the side boundaries, and the distribution of 
grid points along the connecting function. These op- 
tions are for the development of the control functions 

M0> M0> MO. MO. and MO- T yp in s 1 . 2 . 

3, or 4 followed by a CR will cause a display of local 
options relative to each control function. A display 
of local options is shown in figure 11. The first three 
local options define control points within the unit 
square by use of programmed analytical functions 
with the objective of accelerating the process of creat- 
ing relatively simple control functions. Local option 1 
creates a linear distribution of points, local option 2 
creates an exponential distribution of points, and lo- 
cal option 3 creates a distribution from a combination 
of two exponential functions. One user-defined pa- 
rameter is required for local option 2 and two param- 
eters are required for local option 3. Local option 4 
enables the user to digitize points within the unit 
square, thereby defining arbitrary control functions. 
Local option 5 applies only to the bottom and top 
boundaries. It automatically sets the top-boundary 
distribution equal to the bottom-boundary distribu- 
tion, or vice versa. 

Local option 2 generates an exponential distribu- 
tion of points based on the function 

e Kr -1 

t = — ^ (0 < r < 1;0 < f < 1) 

e n — 1 


t = Ki 


e {K 2 /K x )r __ { 
e K 2 — 1 


(0<r < Ai;0<r <Ki) 


r = Ki + (1 + /G) 


e -K 2 (r-K 1 )/(l-K 1 )' 
e~ K * - 1 


(0<Ki< 1 ; 0.001 < K 2 < 50 ) 


Increasingly positive values of If 2 cause decreasing 
first derivative magnitudes of the control function 
near the point (0,0) and the point (1,1) in the unit 
square. An inflection occurs at the point with coordi- 
nates (Ki, K\), where there will be a relatively large 
first derivative magnitude. Conversely, increasingly 
negative values of If 2 cause increasing first derivative 
magnitudes of the control function near points (0,0) 
and (1,1); near the inflection point with coordinates 
(Kx,Ki) the first derivative will be relatively small. 
Figure 13 demonstrates this option. 

After typing local option numbers 2 or 3 followed 
by a CR, the range of the parameter or parameters 
is displayed and the user must respond by typing a 
value or values separated by a comma and followed 
by a CR. 


Arbitrary Control Functions 

As previously stated, local option 4 allows the 
user to specify points within the unit square for con- 
trol function definition. For the five control functions 

MO. MO, MO, M0> and M0> the first P° int is 

always (0,0) and the last point is always (1,1). These 
points are preset by the program. Only the points 
in between are specified by the user, and they must 
be monotonically increasing. The control functions 
fs(0 and /4 (£) have different restrictions and are 
discussed with the orthogonality option (general op- 
tion 5). Typing 4 followed by a CR displays the menu 
and the unit square shown in figure 14. Also, the dig- 
itizing cursor is activated and appears on the screen. 
The digitizing process is the movement of the cursor 
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to place it at desired points inside the unit square 
along with input from the keyboard to signal the 
program to accept a point. In order to implement 
the process, a set of specified commands (keyboard 
characters) must be used, and they are listed on the 
display (fig. 14). The commands are used to inform 
the program on what the user wants done. A de- 
scription of the commands is presented in table III 
to acquaint the user with their purpose. The con- 
trol of the placement of the cursor is accomplished 
through the use of analog devices on the terminal 
(fig. 14). Note that when first digitizing a control 
function, the F function is used to indicate to the 
program the point at abscissa of 1 and the S func- 
tion indicates to the program to proceed to the next 
step. 

Visual Aid for Digitizing Boundary Control 

Functions (P Function) 

The P digitizing function is a visual aid for estab- 
lishing in the unit square points which can be used to 
define the control functions /i(£) and / 2 (C) through 
use of general options 1 and 2. The P function dis- 
plays the boundary curves and draws dotted connect- 
ing grid curves for the number of curves specified at 
the beginning of the interactive session or specified 
with general option 8. It is best to use a small num- 
ber of grid curves to avoid visual confusion. This 
interactive interchange generates points in the unit 
square based on where the user thinks grid curves 
should intersect the boundary curves. 

The option P is invoked by typing the character P 
while digitizing control function /i(£) for the bottom 
boundary and /2(f) for the top boundary (general 
options 1 and 2). This instruction causes the display 
shown in figure 15 to appear on the terminal screen. 
The user is asked which grid line is to be placed at 
a given arc length. The range of grid lines is also 
displayed. The number is typed and followed by a 
CR. The cursor appears on the screen and the user 
places it at the desired point as close as possible to 
the boundary curve. If the grid curve has not been 
used before, the character I is typed to indicate to 
the program that a new point is being digitized. If 
this is a modification to a previously used grid curve, 
the character C is typed to indicate to the program 
that a change to an existing control point is being 
made. Note that a CR does not follow the I or C 
commands. The display of the unit square and asso- 
ciated options appears on the screen with the “loaded 
point” ( 80 i =0, described in eqs. (6)) corresponding 
to what has been digitized. The abscissa of the point 
is I/{N — 1), where I corresponds to the /th grid 
curve and N is the total number of grid curves. The 
ordinate is the normalized approximate arc length 


along the boundary curve nearest the intersection of 
the cross lines. The point is also denoted in the op- 
tion P display for later reference. 

Boundary Distribution Interchange 

A fifth local option is provided for the bottom- 
and top-boundary- grid distributions. This option 
sets the control function for the bottom boundary 
equal to the control function for the top boundary, 
or vice versa. This option is invoked by typing 5 
followed by a CR when it is provided on the menu. 

Control Function Smoothing and Display 

After the control points have been chosen, a dis- 
play of the smooth spline fit to the control points for 
the number of grid points specified at. the beginning 
of the program or in the file RESTART is presented 
on the terminal screen. Three local options relative 
to the spline fit are also presented. An example of 
this display is shown in figure 16. This display al- 
lows the user to ensure that: (1) the control func- 
tion has the desired characteristics; (2) the control 
function, the first derivative, and the second deriva- 
tive are sufficiently continuous; and (3) the control 
function is monotonically increasing for /i(£), /2(C)) 
fo(v)i /6 ( 0 5 an< 3 /V( 0* The spline fit options are: 
(1) change the smoothing parameters tuq and mi in 
equations (6), (2) generate a different control func- 
tion, or (3) accept the distribution and continue to 
display the resultant grid and general options. The 
parameters mo and mi are initially defined in the 
program to be 0.01. 

Grid Orthogonality (General Option 5) 

The two-boundary technique allows control of 
grid orthogonality through the definition of the func- 
tions P(r ) and Q(s) in equations (4), where 

P{r) = K B fs(t) 

Q(s) = K T f 4 {0 

Unlike the previously described grid spacing control 
functions, 7*3 (£) and 74(f) do not have to be mono- 
tonically increasing in the control domain between 
the origin (0,0) and the point (1,1). 

General option 5 is invoked by typing 5 followed 
by a CR when the general options display is on the 
terminal screen. The option causes the display shown 
in figure 17 to appear on the terminal screen. The 
values Kb and K 7- are input from the keyboard as 
alphanumeric constants within the range shown in 
the display. The values are separated by a comma 
and followed by a CR. Three local options are shown 
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in the display. The options are: (0) constant distri- 
butions of fs(C) and /4(£)> (1) user-defined distribu- 
tions of fs(£) and or (2) keep current distri- 

butions. The options are input for both the bottom 
boundary and the top boundary, and the two values 
are separated by a comma and followed by a CR. If lo- 
cal options 0 or 2 are invoked, the grid is recomputed 
and the primary display showing the general options 
and grid appears on the terminal screen. If local op- 
tion 1 is invoked, the display described above for ar- 
bitrary distributions appears on the terminal screen. 
The difference between digitizing control points for 
this option and those previously described is that the 
ordinate of the first point is digitized and the abscissa 
is set equal to 0. When the digitizing option F is ex- 
ecuted, the ordinate of the last point is digitized and 
the abscissa is set equal to 1. The digitizing option S 
indicates to the program to spline smooth the control 
points and to display the smooth control function and 
smoothing options. The digitizing process caused by 
local option 1 is displayed. 

After the control functions fs(£) and f^) have 
been determined, the grid is recomputed and the grid 
and general options are displayed. Note again that 
the orthogonality functions fs(0 and / 4 (£) are ar- 
bitrarily specified according to where the user wants 
the relative effect of the orthogonality to be along the 
bottom and top boundaries. The constants Kq and 
Kx govern the magnitude of the orthogonality. 

Plotting of Grid Derivatives (General 
Option 6) 

An extremely important aspect of an interactive 
grid-generation procedure is the ability to visually as- 
sess the quality of the grid. In addition to inspection 
of the grid itself, visual inspection of the derivatives 
of the physical grid coordinates with respect to the 
computational coordinates is easily made available to 
the interactive user of program TBGG. 

In program TBGG, general option 6 provides for 
the plotting of the scaled grid derivatives along any 
grid curve. General option 6 is invoked by typing 
6 followed by a CR when the general options dis- 
play shown in figure 10 is presented on the terminal 
screen. If this instruction is used, the additional dis- 
play shown in figure 18 also appears in the general 
options display. The additional display asks the user 
to enter the line type (£ or rj coordinate) to be dif- 
ferentiated, the coordinate with respect to which the 
derivative is being taken (1 for x and 2 for ?/), the 
coordinate direction in which the derivative is to be 
plotted (1 for £(/) and 2 for ?/(J)), and the index 
number of the curve along which the derivative is 
taken. The numbers are separated by commas and 
followed by a CR. An example is the derivative of the 


^/-coordinate with respect to the ^-coordinate along 
the 10th curve in the ^-direction. The interactive re- 
sponse is 2,1,2,10 CR. A display similar to the one 
shown at the bottom of figure 18 appears on the ter- 
minal screen. The derivative is plotted normalized 
relative to the lengths of the display axes, and the 
stretch factor is displayed under the plot. The re- 
sponse 0,0, 0,0 CR indicates to the program to end 
the derivative display and causes the general option 
display to appear on the terminal screen. 

Grid Enlargement (General Option 7) 

In order to see the detailed characteristics of a 
grid, it is necessary to plot the grid at a sufficiently 
large scale. On a small screen it is not possible to 
enlarge the entire grid, but a section of the grid can 
be arbitrarily enlarged. General option 7 provides 
for the enlargement of a section of the grid that is 
plotted on the general options display (fig. 10). 

General option 7 is invoked by typing 7 followed 
by a CR. The user is asked the multiple of enlarge- 
ment that is desired. This number is typed followed 
by a CR. The digitizing cross lines appear on the ter- 
minal screen with the present grid. The user places 
the intersection of the cross lines at the center of the 
region that is to be enlarged and types the character 
I. The enlarged grid (fig. 19) appears on the terminal 
screen with the instructions to continue the enlarge- 
ment or to return to the general options. Note that 
a CR does not follow the character I. 

Changing the Number of Grid Points (General 

Option 8) 

A sparse grid can be used in the early stages of 
the development of the control functions and orthog- 
onality parameters that determine the spacing of grid 
points. Using a sparse grid saves time both in the 
computing of the grid and in the drawing of the dis- 
plays, and the amount of information presented to 
the user is not overwhelming. Once the general char- 
acteristics have been determined with use of a sparse 
grid, the number of grid points can be set to any 
desired value. However, to keep the memory size of 
program TBGG within reasonable limits, the present 
maximum grid size is 100 X 100 (100 points in the 
direction and 100 points in the ^-direction). This can 
be changed in the program parameter statements. 
The minimum grid size is 4 x 4. 

General option 8 is invoked by typing 8 followed 
by a CR when the general options and grid (fig. 10) 
are displayed on the terminal screen. A question 
will appear on the terminal screen asking how many 
grid points are desired. The user responds by typ- 
ing the number of desired points in the ^-direction, 
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a comma, the number of points in the //-direction, 
and a CR. The new grid is computed and displayed 
on the terminal screen with the general options. 

Terminating the Program (General Option 9) 

General option 9 terminates the program. It 
also writes the files GRIDOUT and RESTART. Note 
that GRIDOUT and RESTART are local files, and 
the user must actively store them if they are to be 
permanently saved. 

Conclusions 

An algebraic grid-generation technique that is 
easy to understand, easy to apply, and has a high de- 
gree of generality has been presented. Also, an inter- 
active computer program for applying the technique 
in two dimensions has been presented. Boundary def- 
initions are provided as ordered sets of points, and 
the spacing of a grid within the boundaries is gov- 
erned by a set of control functions which are devel- 
oped interactively. The control functions are based 
on an adoption of smooth-cubic-spline functions de- 
fined for a control space, which is the unit square. 


The process allows complete user control of the Her- 
mite interpolation and linear blending. The result is 
an ordered set of points in a physical-coordinate sys- 
tem which corresponds to a uniform grid in a rectan- 
gular computational-coordinate system. Grid deriva- 
tives that are used in the solution of partial differen- 
tial systems can be obtained by numerical differenti- 
ation of the physical grid with respect to the compu- 
tational coordinates and the computation of inverse 
relations. 

The interactive program not only provides a 
means of quickly creating grids, but it also provides 
the opportunity for the user to quickly evaluate and 
draw conclusions about the suitability of grids for 
particular applications. The program is user friendly, 
with prompts for each step of the process, and the 
program tolerates a variety of user errors. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
November 12, 1985 
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Appendix 

Examples 

The first example is the creation of a C-type 
grid about an airfoil. The boundary data for the 
DATANEW file are shown in table IV. The grid, the 
control functions, and the constants Kq and K t are 
shown in figure 20. 

The second example is an O-type grid about 
the same airfoil presented above. In this case the 
bottom-boundary data are the same as those in 
table IV except that they are defined clockwise. The 
“signs” of the orthogonality constants Kg and K t 
are opposite to those in the first example. The top 
boundary is circular and the left and right boundaries 
are identical, extending from the airfoil to the outer 
boundary. Table V shows the DATANEW file for 
this example. Figure 21 shows the grid, the control 
functions, and the constants Kg and Zfj-. 


The third example is a nozzle configuration in 
which no left or right boundaries are specified. The 
boundary data are shown in table VI, and the grid, 
the control functions, and the constants Kg and K j* 
are shown in figure 22. This example is illustrative 
of orthogonality. If the same magnitude of orthogo- 
nality is used in the center of the nozzle as near the 
ends, the grid would overlap. The functions / 3 (C) 
and / 4 (£) both have a value of one at the ends and 
approach zero at the center, which allows the grid to 
have acceptable characteristics. 

The fourth example has the same bottom and top 
boundaries as the third example. Side boundaries 
have been added to show their effect. The functions 
and / 7 (C) control how far into the grid the side 
boundaries modify the initial grid. Table VII shows 
the DATANEW file and figure 23 shows the grid, the 
control functions, and the constants Kq and K y. 
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Symbols and Abbreviations 

A set description of points for smooth- 

ing spline functions 

a, 6, c, d coefficients for smoothed cubic 

splines 


B set of deviations 

C constant specifying extent of 

smoothing 


CR carriage return 


r > derivatives of x- and ^-coordinates 

with respect to r interpolated along 
bottom boundary 


dX j} s ^ , derivatives of x- and ^/-coordinates 

with respect to s interpolated along 
top boundary 


d<t> 


first derivative of cubic-spline 
representation 


d<j>' 2 


second derivative of cubic-spline 
representation 


F 


functional representation of 
smoothed cubic splines 


/l(£), h(£) control functions for distribution of 
grid points along bottom and top 
boundaries 


h{0) A(0 control functions for relative magni- 
tude of normal derivative at bottom 
and top boundaries 

f^(rj) control function for distribution 

of grid points along connecting 
function 


/e(^), /V(£) control functions for relative effect 
of left and right side boundaries 

J, J grid point indices in £ and rj 

directions 


i index for boundary point description 

J Jacobian matrix, 




d£ 1 


’ dx 

dx 

dx ■ 

dx 

dy 

dz 


dS 

dr) 

a* 


dR 




dy 

dy 

dy 

dx 

dy 


a* 

dr) 

a? 

L dx 

§1 

dy 

dz J 


dz 

L ae 

dz. 

dr) 

dz 

aF -1 


3 


index for points along the connect- 
ing function 


K 

Ki,K 2 

K b ,K t 

k 

M 

mo, mi 
N 

Ng, N t 
N l ,N r 
N P 

P(r),Q(s) 
r, s 

f, s 
t 
u 
u 

V , w 
v,vb 

X(I,J), 

X B (r),Y B (r) 


concentration parameter 

control function parameters for bi- 
exponential function 

constants for magnitude of nor- 
mal derivative at bottom and top 
boundaries (KB and KT in com- 
puter displays) 

index denoting point number 
for left- and right-boundary 
descriptions 

total number of points along r\ 
boundary 

parameters to compute deviations 
(MO and Ml in computer displays) 

total number of points along £ 
boundary 

number of points describing bottom 
and top boundaries 

number of points describing left and 
right boundaries 

number of points used to accurately 
compute arc length along connect- 
ing function 

magnitudes of normal derivatives at 
bottom and top boundaries 

normalized approximate arc 
lengths along the bottom and top 
boundaries 

approximate arc lengths along the 
bottom and top boundaries 

parametric variable in Hermite’s 
interpolation functions 

normalized approximate arc length 
along connecting function 

approximate arc length along 
connecting function 

normalized approximate arc lengths 
along left and right boundaries 

approximate arc lengths along left 
and right boundaries 

grid coordinates 

interpolated point along bottom 
boundary 
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Xl,y l 

coordinate description of left 
boundary 

Xr,Yr 

coordinate description of right 
boundary 

X t (s),Y t (s) 

interpolated point along top 
boundary 

x, y,z 

physical-domain coordinates 

a 2 {t), 
a s (t), a 4 (t ) 

interpolation functions for Hermite 
cubic connecting function 

66 

allowed deviation of smoothed 
ordinate from digitized ordinate 

© 1 ) © 2 , 
©3, @4 

parametric set descriptions of left 
and right boundaries 


mO) 

cubic spline representation 


computational-domain coordinates 

T,T 

concentration function variables 

$3, $4 

sets describing bottom and top 
boundaries in parametric form 

4>,e 

abscissa and ordinate of points for 
smoothing splines 


set describing approximate arc 
length along cubic functions 

Notation: 


r\a=c 

\Ja~b 

set of points for enclosed ordered 
pair from initial a value b to final a 
value c 
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TABLE I. DATANEW INPUT DESCRIPTION 


Description 

Card image 

Problem title 

Maximum of 80 alphanumeric 
characters 

Number of grid points* 1 

Bottom boundary, NB=iV^ 

Top boundary, NT=7Vj» 

Left side boundary, NL=iy^ 
Right side boundary, NR =Nr 

NB,NT,NL,NR 

x-coordinates, bottom 
boundary (free format) 

Xi,X2,...,Xnb 

^-coordinates, bottom 
boundary (free format) 

Yi,Y 2 ,...,Y nb 

x-coordinates, top 

boundary (free format) 

X 1 ,X 2 ,...,X N t 

^/-coordinates, top 

boundary (free format) 

Yi,y 2v ..,y nt 

x-coordinates, left 

side boundary (free format) 

Xi,X 2 ,...,Xnl 

^/-coordinates, left 

side boundary (free format) 

Yi,y 2) ...,y nl 

x-coordinates, right 

side boundary (free format) 

X],X 2 ,...,X[mr 

{/-coordinates, right 

side boundary (free format) 

Y 1 ,Y 2 ,...,Y N r 


*The maximum number of points currently allowed for each 
boundary definition is 100. 
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TABLE II. GENERAL OPTIONS 


MB 

Name 

Description 

1 

MODIFY BOTTOM CURVE 

Redefines control function /i(£) for bottom boundary 

2 

MODIFY TOP CURVE 

Redefines control function /2(C) for top boundary 

3 

MODIFY SIDE CURVE 

Redefines control function for concentration effect 

of side boundaries 

4 

MODIFY CONNECTING CURVE 

Redefines control function /s (77) for connecting curve 

5 

CHANGE KB,KT DISTRIBUTION 

Controls orthogonality by assigning values to constants 
Kb and K? and defining control functions fsU) and 

h (0 

6 

PLOT GRID DERIVATIVES 

Displays a derivative along specified grid curve 

7 

BLOW UP A SECTION 

Enlarges grid plot 

8 

CHANGE # OF GRID POINTS 

Assigns number of grid points in each coordinate 
direction 

9 

WRITE OUT GRID AND TERMINATE 

Writes current grid coordinates on file called GRIDOUT, 
creates file RESTART, and ends the session 




TABLE III. DIGITIZING FUNCTIONS 


Option* 

Function 

C 

Corrects a previously digitized point. Place intersection of cross lines near 
point to be corrected and type character C. Move intersection of cross 
lines to desired location and again type character C. New display with 
correction will appear on screen. 

D 

Deletes a point. Place intersection of cross lines near point to be deleted 
and type character D. New display, minus deleted point, will appear 
on screen. 

E 

Erases all points in unit square and allows user to generate completely 
new curve. Type character E. 

F 

Completes definition of a control curve after option E has been invoked 
and new points have been created in unit square. For bottom, top, 
connecting, or side distribution functions, this character indicates to 
program that next and last point is (1,1). For orthogonality options, 
this character indicates that next and last point is (1 ,y), where y is 
y-coordinate of intersection of cross lines. Type character F. 

I 

Digitizes point in unit square where cross lines intersect. Type charac- 
ter I. 

L 

Forces spline function to pass exactly through a point. Place intersection 
of cross lines at point to be digitized and type character L. Equa- 
tions (6) are bypassed and deviation 60(1) is set to zero. Such points 
are referred to as “loaded points.” Type character L. 

P 

Replaces unit square display with a display of boundary curves. This 
additional display allows input of a grid point number and a digitized 
position on a boundary curve to generate control points in unit square. 
This option only applies to bottom and top boundaries. Type charac- 
ter P. 

S 

Indicates to program that digitizing process is complete or there are no 
changes. This option indicates to program to go to next step. Type 
character S. 


These commands are not followed by a CR. 




TABLE IV. DATANEW FILE FOR AIRFOIL C-TYPE GRID EXAMPLE 


AIRFOIL C-TYPE GRID EXAMPLE 
81 39 2 2 


1 .009020 

1.007465 

1.002809 

.995081 

.934676 

.912667 

.888143 

.861253 

.697578 

.660412 

.622286 

.583433 

.386735 

.348608 

.311442 

.275467 

. 120878 

.096353 

.074345 

. 054988 

.001555 

. 000000 

.001555 

.006211 

.074345 

.096353 

. 120878 

.147768 

.311442 

.348608 

.386735 

.425587 

.622286 

.660412 

.697578 

.733553 

.888143 

.912667 

.934676 

.954032 

1 .007465 

1 . 009020 



.000000 

.000208 

.000865 

.001948 

.010071 

.012888 

.015943 

.019195 

.036874 

.040379 

.043762 

.046972 

.058497 

.059514 

.059987 

.059878 

. 050002 

.046201 

.041854 

.037006 

.006907 

.000000 

-.006907 

-.013562 

- .041854 

-.046201 

-.050002 

- . 053222 

-.059987 

-.059514 

-.058497 

-.056980 

- .043762 

-.040379 

-.036874 

-.033297 

-.015943 

-.012888 

-.010071 

-.007532 

- .000208 

.000000 



1 .00902 

.000000 

-.087156 

-.173648 

-.573576 

-.642788 

- .707107 

-.766044 

- .965926 

- . 984808 

-.996195 - 

1.000000 

-.906308 

-.866025 

-.819152 

-.766044 

- .422618 

- . 342020 

-.258819 

-.173648 

1.000000 

1.000000 

.996195 

.984808 

.819152 

.766044 

.707107 

.642788 

.258819 

.173648 

.087156 

.000000 

- .422618 

- . 500000 

-.573576 

-.642788 

-.906308 

-.939693 

-.965926 

-.984808 

1 .00902 

1.00902 




0.0 1.0 

1.00902 1.00902 

0.0 - 1.0 


.984328 

.970617 

.954032 


.832163 

.801054 

.768116 

.733553 

.544094 

.504510 

.464927 

.425587 

.240904 

.207967 

.176857 

.147768 

.038404 

.024692 

.013940 

.006211 

.013940 

.024692 

.038404 

.054988 

.176857 

.207967 

.240904 

.275467 

.464927 

.504510 

.544094 

.583433 

.768116 

.801054 

.832163 

.861253 

.970617 

.984328 

.995081 

1.002809 

.003438 

.005309 

.007532 


.022600 

.026114 

.029695 

.033297 

.049954 

.052653 

.055014 

.056980 

.059159 

.057812 

.055831 

.053222 

.031705 

.026000 

.019938 

.013562 

-.019938 

-.026000 

-.031705 

-.037006 

-.055831 

-.057812 

- .059159 

- .059878 

-.055014 

- . 052653 

-.049954 

-.046972 

-.029695 

-.026114 

-.022600 

-.019195 

-.005309 

-.003438 

-.001948 

-.000865 

- .258819 

-.342020 

-.422618 

-.500000 

-.819152 

-.866025 

-.906308 

-.939693 

-.996195 

-.984808 

- .965926 

- . 939693 

-.707107 

-.642788 

-.573576 

-.500000 

-.087156 

.000000 

1.00902 


.965926 

.939693 

.906308 

.866025 

.573576 

.500000 

.422618 

.342020 

-.087156 

-.173648 

-.258819 

-.342020 

-.707107 

-.766044 

-.819152 

-.866025 

-.996195 

- 1.000000 

- 1.00000 
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TABLE V. DATANEW FILE FOR AIRFOIL O-TYPE GRID EXAMPLE 


AIRFOIL O-TYPE GRID EXAMPLE 
81 91 2 2 


1.009020 

1.007465 

1.002809 

.995081 

.984328 

.970617 

.954032 

.934676 

.912667 

.888143 

.861253 

.832163 

.801054 

.768116 

.733553 

.697578 

.660412 

.622286 

.583433 

.544094 

.504510 

.464927 

.425587 

.386735 

.348608 

.311442 

.275467 

.240904 

.207967 

.176857 

.147768 

.120878 

.096353 

.074345 

.054988 

.038404 

.024692 

.013940 

.006211 

.001555 

.000000 

.001555 

.006211 

.013940 

.024692 

.038404 

.054988 

.074345 

.096353 

.120878 

. 147768 

. 176857 

.207967 

.240904 

.275467 

.311442 

.348608 

.386735 

.425587 

.464927 

.504510 

.544094 

.583433 

.622286 

.660412 

.697578 

.733553 

.768116 

.801054 

.832163 

.861253 

.888143 

.912667 
1 .009020 

.934676 

.954032 

.970617 

.984328 

.995081 

1.002809 

1.007465 

.000000 

- .000208 

- .000865 

- .001948 

- .003438 

- .005309 

- .007532 

-.010071 

- . 012888 

- .015943 

- .019195 

-.022600 

- .026114 

- .029695 

-.033297 

-.036874 

- .040379 

- .043762 

- .046972 

- .049954 

- .052653 

-.055014 

-.056980 

- .058497 

- .059514 

- .059987 

- .059878 

-.059159 

-.057812 

- .055831 

-.053222 

- .050002 

- .046201 

-.041854 

- .037006 

- .031705 

-.026000 

- .019938 

-.013562 

- .006907 

.000000 

.006907 

.013562 

.019938 

.026000 

.031705 

.037006 

.041854 

.046201 

.050002 

.053222 

.055831 

.057812 

.059159 

.059878 

.059987 

.059514 

.058497 

.056980 

.055014 

.052653 

.049954 

.046972 

.043762 

.040379 

.036874 

.033297 

.029695 

.026114 

.022600 

.019195 

.015943 

.012888 
. oooooo 

.010071 

.007532 

.005309 

.003438 

.001948 

.000865 

.000208 

2.000000 

1.996346 

1.985402 

1.967221 

1.941893 

1.909539 

1.870318 

1.824421 

1.772072 

1.713525 

1 . 649067 

1.579010 

1.503696 

1.423492 

1.338789 

1.250000 

1.157557 

1.061910 

.963525 

. 862883 

.760472 

.656793 

.552349 

.447651 

.343207 

.239528 

.137117 

.036475 

-.061910 

-.157557 

-.250000 

- .338789 

- .423492 

-.503696 

-.579010 

- . 649067 

-.713525 

- . 772072 

-.824421 

-.870318 

- .909539 

-.941893 

-.967221 

-.985402 

- .996346 

-1.000000 

-.996346 

-.985402 

- .967221 

- .941893 

- .909539 

- .870318 

- .824421 

- .772072 

-.713525 

-.649067 

- . 579010 

- .503696 

-.423492 

-.338789 

-.250000 

-.157557 

-.061910 

.036475 

.137117 

.239528 

.343207 

.447651 

.552349 

.656793 

.760472 

.862883 

.963525 

1.061910 

1 . 157557 

1.250000 

1.338789 

1.423492 

1.503696 

1.579010 

1.649067 

1.985402 

1.713525 

1.996346 

1.772072 

2.000000 

1.824421 

1.870318 

1.909539 

1.941893 

1.967221 

.000000 

- . 104635 

-.208760 

-.311868 

- .413456 

-.513030 

-.610105 

- .704207 

-.794879 

-.881678 

-.964181 - 

1.041988 - 

1.114717 

-1.182016 - 

1.243556 - 

1.299038 

-1.348191 

-1.390776 

-1.426585 

-1.455444 

-1.477212 

-1.491783 

-1.499086 

-1.49908 

-1.491783 

-1.477212 

-1.455444 

-1.426585 

-1.390776 

-1.348191 

-1.299038 

-1.24355 

-1.182016 

-1.114717 

-1.041988 

- .964181 

-.881678 

-.794879 

-.704207 

-.61010 

- .513030 

-.413456 

-.311868 

-.208760 

-.104635 

.OOOOOO 

.104635 

.208760 

.311868 

.413456 

.513030 

.610105 

.704207 

.794879 

.881678 

.964181 

1.041988 

1.114717 

1.182016 

1.243556 

1.299038 

1.348191 

1.390776 

1.426585 

1 . 455444 

1.477212 

1.491783 

1.499086 

1.499086 

1.491783 

1.477212 

1 .455444 

1.426585 

1.390776 

1.348191 

1.299038 

1 . 243556 

1.182016 

1 .114717 

1 .041988 

.964181 
.208760 
1 .009020 
.000000 
1.009020 
. OOOOOO 

.881678 
. 104635 
2.000000 
.OOOOOO 
2.000000 
.OOOOOO 

.794879 

.OOOOOO 

.704207 

.610105 

.513030 

.413456 

.311868 
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TABLE VI. DATANEW FILE FOR NOZZLE CONFIGURATION MINUS 


SIDE BOUNDARIES 

NOZZLE CONFIGURATION MINUS SIDE BOUNDARIES 
50 50 0 0 


- 13.00000 

- 10.34694 

- 7.69388 

- 5.04082 

- 2.38776 

.26531 

2.91837 

5.57143 

8.22449 

10.87755 

- 7.00000 

- 5.70148 

- 4.44547 

- 3.28115 

- 2.34907 

- 2.00468 

- 2.50356 

- 3.50219 

- 4.69161 

- 5.95866 

- 13.00000 

- 10.34694 

- 7.69388 

- 5.04082 

- 2.38776 

.26531 

2.91837 

5.57143 

8.22449 

10.87755 

7.00000 

5.70148 

4.44547 

3.28115 

2.34907 

2.00468 

2.50356 

3.50219 

4.69161 

5.95866 


- 12.46939 
- 9.81633 
- 7.16327 
- 4.51020 
- 1.85714 
.79592 
3.44898 
6.10204 
8.75510 
11.40816 
- 6.73806 
- 5.44592 
- 4.20274 
- 3.06863 
- 2.21774 
- 2.04173 
- 2.67721 
- 3.73023 
- 4.94067 
- 6.21726 
- 12.46939 
- 9.81633 
- 7.16327 
- 4.51020 
- 1.85714 
.79592 
3.44898 
6.10204 
8.75510 
11.40816 
6.73806 
5.44592 
4.20274 
3.06863 
2.21774 
2.04173 
2.67721 
3.73023 
4 . 94067 
6.21726 


- 11.93878 

- 9.28571 

- 6.63265 

- 3.97959 

- 1.32653 

1.32653 

3.97959 

6.63265 

9.28571 

11.93878 

- 6.47711 

- 5.19222 

- 3.96407 

- 2.86653 

- 2.11390 

- 2.11390 

- 2.86653 

- 3.96407 

- 5.19222 

- 6.47711 

- 11.93878 

- 9.28571 

- 6.63265 

- 3.97959 

- 1.32653 

1.32653 

3.97959 

6.63265 

9.28571 

11.93878 

6.47711 

5.19222 

3.96407 

2.86653 

2.11390 

2.11390 

2.86653 

3.96407 

5.19222 

6.47711 


- 11.40816 
- 8.75510 
- 6.10204 
- 3.44898 
-.79592 
1.85714 
4.51020 
7.16327 
9.81633 
12.46939 
- 6.21726 
- 4.94067 
- 3.73023 
- 2.67721 
- 2.04173 
- 2.21774 
- 3.06863 
- 4.20274 
- 5.44592 
- 6.73806 
- 11.40816 
- 8.75510 
- 6.10204 
- 3.44898 
- .79592 
1.85714 
4.51020 
7.16327 
9.81633 
12.46939 
6.21726 
4.94067 
3.73023 
2.67721 
2.04173 
2.21774 
3.06863 
4.20274 
5.44592 
6.73806 


- 10.87755 
- 8.22449 
- 5.57143 
- 2.91837 
-.26531 
2.38776 
5.04082 
7.69388 
10.34694 

13.00000 
- 5.95866 
- 4.69161 
- 3.50219 
- 2.50356 
- 2.00468 
- 2.34907 
- 3.28115 
- 4.44547 
- 5.70148 
- 7.00000 
- 10.87755 
- 8.22449 
- 5.57143 
- 2.91837 
- .26531 
2.38776 
5.04082 
7.69388 
10.34694 

13.00000 
5.95866 
4.69161 
3.50219 
2.50356 
2.00468 
2.34907 
3.28115 
4.44547 
5.70148 

7,00000 



TABLE VII. DATANEW FILE FOR NOZZLE CONFIGURATION 
WITH SIDE BOUNDARIES 


NOZZLE CONFIGURATION 

WITH SIDE 

BOUNDARIES 


50 50 21 

21 




-13.00000 

-12.46939 - 

-11.93878 

-11.40816 

-10.87755 

-10.34694 

-9.81633 

-9.28571 

-8.75510 

-8.22449 

-7.69388 

-7.16327 

-6.63265 

-6.10204 

-5.57143 

-5.04082 

-4.51020 

-3.97959 

-3.44898 

-2.91837 

-2.38776 

-1.85714 

-1.32653 

-.79592 

-.26531 

.26531 

.79592 

1.32653 

1.85714 

2.38776 

2.91837 

3.44898 

3.97959 

4.51020 

5.04082 

5.57143 

6.10204 

6.63265 

7.16327 

7.69388 

8.22449 

8.75510 

9.28571 

9.81633 

10.34694 

10.87755 

11.40816 

11.93878 

12.46939 

13.00000 

-7.00000 

-6.73806 

-6.47711 

-6.21726 

-5.95866 

-5.70148 

-5.44592 

-5.19222 

-4.94067 

-4.69161 

-4.44547 

-4.20274 

-3.96407 

-3.73023 

-3.50219 

-3.28115 

-3.06863 

-2.86653 

-2.67721 

-2.50356 

-2.34907 

-2.21774 

-2.11390 

-2.04173 

-2.00468 

-2.00468 

-2.04173 

-2.11390 

-2.21774 

-2.34907 

-2.50356 

-2.67721 

-2.86653 

-3.06863 

-3.28115 

-3.50219 

-3.73023 

-3.96407 

-4.20274 

-4.44547 

-4.69161 

-4.94067 

-5.19222 

-5.44592 

-5.70148 

-5.95866 

-6.21726 

-6.47711 

-6.73806 

-7.00000 

-13.00000 

-12.46939 - 

-11.93878 

-11.40816 

-10.87755 

-10.34694 

-9.81633 

-9.28571 

-8.75510 

-8.22449 

-7.69388 

-7.16327 

-6.63265 

-6.10204 

-5.57143 

-5.04082 

-4.51020 

-3.97959 

-3.44898 

-2.91837 

-2.38776 

-1.85714 

-1.32653 

-.79592 

-.26531 

.26531 

.79592 

1.32653 

1.85714 

2.38776 

2.91837 

3.44898 

3 . 97959 

4.51020 

5.04082 

5.57143 

6.10204 

6.63265 

7.16327 

7.69388 

8.22449 

8.75510 

9.28571 

9.81633 

10.34694 

10.87755 

11.40816 

11.93878 

12.46939 

13.00000 

7.00000 

6.73806 

6.47711 

6.21726 

5.95866 

5.70148 

5.44592 

5.19222 

4.94067 

4.69161 

4.44547 

4.20274 

3.96407 

3.73023 

3.50219 

3.28115 

3.06863 

2 . 86653 

2.67721 

2.50356 

2.34907 

2.21774 

2.11390 

2.04173 

2.00468 

2.00468 

2.04173 

2.11390 

2.21774 

2.34907 

2.50356 

2.67721 

2.86653 

3.06863 

3.28115 

3.50219 

3.73023 

3.96407 

4.20274 

4.44547 

4.69161 

4.94067 

5.19222 

5.44592 

5.70148 

5.95866 

6.21726 

6.47711 

6.73806 

7.00000 

13.0 -12.073 -11. 2366 - 

10.573 -10 

.1467 


10.0 -10.1467 -10.573 - 

11.2366 -12.073 -13.0 


13.927 -14. 

7625 -15.427 

' -15.8533 

-16.0 


15.8533 -15 

.427 -14.7625 -13.927 

-13.0 


7.0 -6.29999 -5.59999 - 

4.89999 -4 

.19999 


3.49999 -2. 

80001 -2.10001 -1.39996 -0.7001 


.0 0.7001 1 

.39996 2.10001 2.80001 

3.49999 


.19999 4.69999 5.59999 

6.29999 7.' 

0 


13.00000 

12.45946 

11.95849 

11.50226 

11.09628 

10.74625 

10.45780 

10.23613 

10.08564 

10.00955 

10.00955 

10.08564 

10.23613 

10.45780 

10.74625 

11.09628 

11.50226 

11.95849 

12.45946 

13.00000 

-7.00000 

-6.26316 

-5.52632 

-4.78947 

-4.05263 

-3.31579 

-2.57895 

-1.84211 

-1.10526 

- .36842 

.36842 

1 .10526 

1.84211 

2.57895 

3.31579 

4.05263 

4.78947 

5.52632 

6.26316 

7.00000 







(0,0) (1,0) 


Computational coordinate 
Figure 4. Control domain. 
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Figure 5. Cubic blending functions. 


Control functions 



Function evaluation 
plus 

interpolation 

















Figure 9. Smooth-spline control functions. 


OPTIONS 

1. MODIFY BOTTOM CURVE 

2. MODIFY TOP CURVE 

3. MODIFY SIDE CURVE 

4. MODIFY CONNECTING CURVE 

5. CHANGE KB.KT DISTRIBUTION 
ENTER INTEGER VALUE AND PRESS 


6. PLOT GRID DERIVATIVES 

7. BLOW UP A SECTION 

8. CHANGE # OF GRID POINTS 

9. WRITE OUT GRID AND TERMINATE 


"RETURN" 



Figure 10. General options and grid display. 
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OPTIONS 

ENTER INTEGER VALUE FOR TYPE OF NOMINAL CURVE 

1 LINEAR 

2 SINGLE EXPONENTIAL 

3 BI-EXPONENTIAL 

4 USER DEFINED 

5 EQUATE TO TOP CURVE (FOR BOTTOM CURVE) 


Figure 11. Local option display for grid spacing control. 


( 0 , 1 ) ( 1 , 



Figure 12. Single exponential function for grid spacing control. 






OPTIONS FUNCTIONS 

C CORRECT A PREVIOUSLY DIGITIZED POINT 

D DELETE A POINT 

E ERASE ALL POINTS 

F COMPLETE THE DEFINITION OF A CONTROL CURVE 
I DIGITIZE A POINT IN THE UNIT SQUARE 
L LOAD A POINT (NO SMOOTHING) 

P PRESENT BOUNDARY CURVE DISPLAY 
S COMPLETE DIGITIZING PROCESS 
TYPE LETTER FOR DESIRED FUNCTION 



Figure 14. Display for creating arbitrary control functions. 



THE BOUNDARY BEING WORKED ON IS A SOLID LINE 
GRID LINE #1 IS MARKED BY A DIAMOND 
LOADED POINTS ARE SHOWN BY A “+" SIGN 
A MAXIMUM OF 10 POINTS CAN BE LOADED 
CURRENTLY 1 POINT IS LOADED 


THE LOAD OPTIONS ARE: 

1: SHIFT AN ORIGINAL GRID LINE TO A NEW LOCATION. RANGE 2 THRU 10 
2: CHANGE A SHIFTED GRID LINE POSITION 
ENTER THE GRID LINE NUMBER (ENTER 0 TO STOP) AND PRESS ' RETURN" 


LETTER FUNCTION 

C TO CHANGE A SHIFTED GRID LINE 

I TO SHIFT A GRID LINE 

S TO STOP LOADING MODE 


MOVE CURSOR TO 
POSITION ON BOUNDARY 
THEN TYPE LETTER 
FOR DESIRED FUNCTION 



Figure 15. Display to aid in control of grid spacing along boundaries. 



OPTIONS 

1 CHANGE THE SMOOTHING OF THE CURVE 

2 CHANGE THE CURVE 

3 NO CHANGE 

ENTER INTEGER VALUE AND PRESS "RETURN" 

"1" CURVE SMOOTHING OPTION 

THE VARIABLES MO AND Ml ENABLE THE USER TO SMOOTH THE CURVE 

GENERATED BY THE PROGRAM 

THE CURRENT VALUES OF MO AND Ml ARE .010. .010 

TYPE IN THE VALUES SEPARATED BY A COMMA 



Function 
2nd derivative 

1st derivative 


Figure 16. Display of smooth-grid-spacing control function and local options. 


SELECT MAGNITUDES OF ORTHOGONALITY OF GRID 
KB (ABSOLUTE RANGE 0.0 TO 30.0) 

KT (ABSOLUTE RANGE 0.0 TO 15.0) 

PRESENT VALUES ARE .000. .000 

ENTER TWO REAL NUMBERS SEPARATED BY A COMMA AND PRESS "RETURN" 

SELECT KB.KT DISTRIBUTIONS; OPTIONS ARE: 

0 CONSTANT DISTRIBUTION 

1 USER-DEFINED DISTRIBUTION 

2 KEEP CURRENT DISTRIBUTION 
CURRENT VALUES ARE 0.0 

ENTER TWO INTEGER VALUES SEPARATED BY A COMMA AND PRESS "RETURN" 

Figure 17. Display for grid orthogonality (general option 5). 



SELECT LINE TYPE, LINE NUMBER, DERIVATIVE TYPE, AND DIRECTION 
LINE TYPE 1 ALONG BOUNDARY 

2 ACROSS BOUNDARIES 
DERIVATIVE TYPE 1 FOR X 
2 FOR Y 

DERIVATIVE DIRECTION 1 ALONG BOUNDARY 

2 ACROSS BOUNDARIES 
ENTER 0,0, 0,0 TO QUIT THIS MODE 

ENTER FOUR INTEGER VALUES SEPARATED BY COMMAS AND PRESS "RETURN " 



Figure 18. Display for grid derivatives. 



SELECT THE MULTIPLE TO CHANGE THE SIZE OF THE GRID BY 
(1.0 TO RESTORE THE ORIGINAL GRID. 0.0 TO EXIT MODE) 
ENTER REAL VALUE AND PRESS "RETURN" 


PLACE CURSOR AT THE CENTER OF THE REGION TO BE BLOWN UP AND PRESS 
THE LETTER T 



Figure 19. Display for grid enlargement (general option 7). 
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